library(BIEN)
Loading required package: RPostgreSQL
Loading required package: DBI
Type vignette("BIEN") to get started
The BIEN database has been updated to version 4.1 (as of 31-10-2018)!  This version includes nearly twice as many occurrence records and improved validations.
library(maps) 
library(sp)
library(raster) 
vignette("BIEN")
A_racemosa <- BIEN_occurrence_species(species = "Actaea racemosa")
M_canadense <- BIEN_occurrence_species(species = "Menispermum canadense")
V_prunifolium <- BIEN_occurrence_species(species = "Viburnum prunifolium")

A_racemosa_rangemap <- BIEN_ranges_load_species("Actaea racemosa")
M_canadense_rangemap <- BIEN_ranges_load_species("Menispermum canadense")
V_prunifolium_rangemap <- BIEN_ranges_load_species("Viburnum prunifolium")
# Both can be mapped, for example: 
map('world', fill=T, col="grey", bg="light blue", xlim=c(-145,-50), ylim=c(12,78))
  plot(M_canadense_rangemap, col="green", add=T)

# Map another Jungia species
map('world', fill=T, col="grey", bg="light blue", xlim=c(-145,-50),ylim=c(12,78))
plot(V_prunifolium_rangemap, col="red", add=T)

# And a third Jungia
map('world', fill=T, col="grey", bg="light blue", xlim=c(-145,-50),ylim=c(12,78))
plot(A_racemosa_rangemap, col="green", add=T)

# Map occurence records for three species of Jungia
map('world', fill=T, col="grey", bg="light blue", xlim=c(-145,-50),ylim=c(12,78))
points(cbind(M_canadense$longitude, M_canadense$latitude), col="red", pch=20, cex=.8)
points(cbind(A_racemosa$longitude, A_racemosa$latitude), col="blue", pch=20, cex=.8)
points(cbind(V_prunifolium$longitude, V_prunifolium$latitude), col="green", pch=20, cex=.8)

all.worldclim <- raster::getData("worldclim", res = 10, var = "bio")
trying URL 'https://biogeo.ucdavis.edu/data/climate/worldclim/1_4/grid/cur/bio_10m_bil.zip'
Content type 'application/zip' length 10735619 bytes (10.2 MB)
==================================================
downloaded 10.2 MB
all.worldalt <- raster::getData("worldclim", res = 10, var = "alt")
trying URL 'https://biogeo.ucdavis.edu/data/climate/worldclim/1_4/grid/cur/alt_10m_bil.zip'
Content type 'application/zip' length 933921 bytes (912 KB)
==================================================
downloaded 912 KB
plot(all.worldclim[["bio1"]])

plot(all.worldclim[["bio12"]])

plot(all.worldalt[["alt"]])

plot(all.worldclim[["bio1"]])
points(V_prunifolium[,c("longitude", "latitude")], pch = 24, cex = 0.5)

plot(all.worldclim[["bio12"]])
points(V_prunifolium[,c("longitude", "latitude")], pch = 16, cex = 0.5, col = "green")

plot(all.worldalt[["alt"]], xlim=c(-145,-50),ylim=c(12,78))
points(V_prunifolium[,c("longitude", "latitude")], pch = 16, cex = 0.5, col = "green")
points(A_racemosa[,c("longitude", "latitude")], pch = 16, cex = 0.5, col = "blue")
points(M_canadense[,c("longitude", "latitude")], pch = 16, cex = 0.5, col = "red")

plot(all.worldclim[["bio1"]])
points(A_racemosa[,c("longitude", "latitude")], pch = 24, cex = 0.5)

plot(all.worldclim[["bio12"]])
points(A_racemosa[,c("longitude", "latitude")], pch = 16, cex = 0.5, col = "blue")

plot(all.worldalt[["alt"]], xlim=c(-145,-50),ylim=c(12,78))
points(A_racemosa[,c("longitude", "latitude")], pch = 16, cex = 0.5, col = "blue")

plot(all.worldclim[["bio1"]])
points(M_canadense[,c("longitude", "latitude")], pch = 24, cex = 0.5)

plot(all.worldclim[["bio12"]])
points(M_canadense[,c("longitude", "latitude")], pch = 16, cex = 0.5, col = "red")

plot(all.worldalt[["alt"]], xlim=c(-145,-50),ylim=c(12,78))
points(M_canadense[,c("longitude", "latitude")], pch = 16, cex = 0.5, col = "red")

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CmxpYnJhcnkoQklFTikKbGlicmFyeShtYXBzKSAKbGlicmFyeShzcCkKbGlicmFyeShyYXN0ZXIpIAp2aWduZXR0ZSgiQklFTiIpCmBgYAoKYGBge3J9CkFfcmFjZW1vc2EgPC0gQklFTl9vY2N1cnJlbmNlX3NwZWNpZXMoc3BlY2llcyA9ICJBY3RhZWEgcmFjZW1vc2EiKQpNX2NhbmFkZW5zZSA8LSBCSUVOX29jY3VycmVuY2Vfc3BlY2llcyhzcGVjaWVzID0gIk1lbmlzcGVybXVtIGNhbmFkZW5zZSIpClZfcHJ1bmlmb2xpdW0gPC0gQklFTl9vY2N1cnJlbmNlX3NwZWNpZXMoc3BlY2llcyA9ICJWaWJ1cm51bSBwcnVuaWZvbGl1bSIpCgpBX3JhY2Vtb3NhX3JhbmdlbWFwIDwtIEJJRU5fcmFuZ2VzX2xvYWRfc3BlY2llcygiQWN0YWVhIHJhY2Vtb3NhIikKTV9jYW5hZGVuc2VfcmFuZ2VtYXAgPC0gQklFTl9yYW5nZXNfbG9hZF9zcGVjaWVzKCJNZW5pc3Blcm11bSBjYW5hZGVuc2UiKQpWX3BydW5pZm9saXVtX3JhbmdlbWFwIDwtIEJJRU5fcmFuZ2VzX2xvYWRfc3BlY2llcygiVmlidXJudW0gcHJ1bmlmb2xpdW0iKQpgYGAKCmBgYHtyfQojIEJvdGggY2FuIGJlIG1hcHBlZCwgZm9yIGV4YW1wbGU6IAptYXAoJ3dvcmxkJywgZmlsbD1ULCBjb2w9ImdyZXkiLCBiZz0ibGlnaHQgYmx1ZSIsIHhsaW09YygtMTQ1LC01MCksIHlsaW09YygxMiw3OCkpCiAgcGxvdChNX2NhbmFkZW5zZV9yYW5nZW1hcCwgY29sPSJncmVlbiIsIGFkZD1UKQpgYGAKCmBgYHtyfQojIE1hcCBhbm90aGVyIEp1bmdpYSBzcGVjaWVzCm1hcCgnd29ybGQnLCBmaWxsPVQsIGNvbD0iZ3JleSIsIGJnPSJsaWdodCBibHVlIiwgeGxpbT1jKC0xNDUsLTUwKSx5bGltPWMoMTIsNzgpKQpwbG90KFZfcHJ1bmlmb2xpdW1fcmFuZ2VtYXAsIGNvbD0icmVkIiwgYWRkPVQpCmBgYAoKYGBge3J9CiMgQW5kIGEgdGhpcmQgSnVuZ2lhCm1hcCgnd29ybGQnLCBmaWxsPVQsIGNvbD0iZ3JleSIsIGJnPSJsaWdodCBibHVlIiwgeGxpbT1jKC0xNDUsLTUwKSx5bGltPWMoMTIsNzgpKQpwbG90KEFfcmFjZW1vc2FfcmFuZ2VtYXAsIGNvbD0iZ3JlZW4iLCBhZGQ9VCkKYGBgCgpgYGB7cn0KIyBNYXAgb2NjdXJlbmNlIHJlY29yZHMgZm9yIHRocmVlIHNwZWNpZXMgb2YgSnVuZ2lhCm1hcCgnd29ybGQnLCBmaWxsPVQsIGNvbD0iZ3JleSIsIGJnPSJsaWdodCBibHVlIiwgeGxpbT1jKC0xNDUsLTUwKSx5bGltPWMoMTIsNzgpKQpwb2ludHMoY2JpbmQoTV9jYW5hZGVuc2UkbG9uZ2l0dWRlLCBNX2NhbmFkZW5zZSRsYXRpdHVkZSksIGNvbD0icmVkIiwgcGNoPTIwLCBjZXg9LjgpCnBvaW50cyhjYmluZChBX3JhY2Vtb3NhJGxvbmdpdHVkZSwgQV9yYWNlbW9zYSRsYXRpdHVkZSksIGNvbD0iYmx1ZSIsIHBjaD0yMCwgY2V4PS44KQpwb2ludHMoY2JpbmQoVl9wcnVuaWZvbGl1bSRsb25naXR1ZGUsIFZfcHJ1bmlmb2xpdW0kbGF0aXR1ZGUpLCBjb2w9ImdyZWVuIiwgcGNoPTIwLCBjZXg9LjgpCmBgYAoKYGBge3J9CmFsbC53b3JsZGNsaW0gPC0gcmFzdGVyOjpnZXREYXRhKCJ3b3JsZGNsaW0iLCByZXMgPSAxMCwgdmFyID0gImJpbyIpCmFsbC53b3JsZGFsdCA8LSByYXN0ZXI6OmdldERhdGEoIndvcmxkY2xpbSIsIHJlcyA9IDEwLCB2YXIgPSAiYWx0IikKcGxvdChhbGwud29ybGRjbGltW1siYmlvMSJdXSkKcGxvdChhbGwud29ybGRjbGltW1siYmlvMTIiXV0pCnBsb3QoYWxsLndvcmxkYWx0W1siYWx0Il1dKQpgYGAKCmBgYHtyfQpwbG90KGFsbC53b3JsZGNsaW1bWyJiaW8xIl1dKQpwb2ludHMoVl9wcnVuaWZvbGl1bVssYygibG9uZ2l0dWRlIiwgImxhdGl0dWRlIildLCBwY2ggPSAyNCwgY2V4ID0gMC41KQpwbG90KGFsbC53b3JsZGNsaW1bWyJiaW8xMiJdXSkKcG9pbnRzKFZfcHJ1bmlmb2xpdW1bLGMoImxvbmdpdHVkZSIsICJsYXRpdHVkZSIpXSwgcGNoID0gMTYsIGNleCA9IDAuNSwgY29sID0gImdyZWVuIikKcGxvdChhbGwud29ybGRhbHRbWyJhbHQiXV0sIHhsaW09YygtMTQ1LC01MCkseWxpbT1jKDEyLDc4KSkKcG9pbnRzKFZfcHJ1bmlmb2xpdW1bLGMoImxvbmdpdHVkZSIsICJsYXRpdHVkZSIpXSwgcGNoID0gMTYsIGNleCA9IDAuNSwgY29sID0gImdyZWVuIikKcG9pbnRzKEFfcmFjZW1vc2FbLGMoImxvbmdpdHVkZSIsICJsYXRpdHVkZSIpXSwgcGNoID0gMTYsIGNleCA9IDAuNSwgY29sID0gImJsdWUiKQpwb2ludHMoTV9jYW5hZGVuc2VbLGMoImxvbmdpdHVkZSIsICJsYXRpdHVkZSIpXSwgcGNoID0gMTYsIGNleCA9IDAuNSwgY29sID0gInJlZCIpCmBgYApgYGB7cn0KcGxvdChhbGwud29ybGRjbGltW1siYmlvMSJdXSkKcG9pbnRzKEFfcmFjZW1vc2FbLGMoImxvbmdpdHVkZSIsICJsYXRpdHVkZSIpXSwgcGNoID0gMjQsIGNleCA9IDAuNSkKcGxvdChhbGwud29ybGRjbGltW1siYmlvMTIiXV0pCnBvaW50cyhBX3JhY2Vtb3NhWyxjKCJsb25naXR1ZGUiLCAibGF0aXR1ZGUiKV0sIHBjaCA9IDE2LCBjZXggPSAwLjUsIGNvbCA9ICJibHVlIikKcGxvdChhbGwud29ybGRhbHRbWyJhbHQiXV0sIHhsaW09YygtMTQ1LC01MCkseWxpbT1jKDEyLDc4KSkKcG9pbnRzKEFfcmFjZW1vc2FbLGMoImxvbmdpdHVkZSIsICJsYXRpdHVkZSIpXSwgcGNoID0gMTYsIGNleCA9IDAuNSwgY29sID0gImJsdWUiKQpgYGAKYGBge3J9CnBsb3QoYWxsLndvcmxkY2xpbVtbImJpbzEiXV0pCnBvaW50cyhNX2NhbmFkZW5zZVssYygibG9uZ2l0dWRlIiwgImxhdGl0dWRlIildLCBwY2ggPSAyNCwgY2V4ID0gMC41KQpwbG90KGFsbC53b3JsZGNsaW1bWyJiaW8xMiJdXSkKcG9pbnRzKE1fY2FuYWRlbnNlWyxjKCJsb25naXR1ZGUiLCAibGF0aXR1ZGUiKV0sIHBjaCA9IDE2LCBjZXggPSAwLjUsIGNvbCA9ICJyZWQiKQpwbG90KGFsbC53b3JsZGFsdFtbImFsdCJdXSwgeGxpbT1jKC0xNDUsLTUwKSx5bGltPWMoMTIsNzgpKQpwb2ludHMoTV9jYW5hZGVuc2VbLGMoImxvbmdpdHVkZSIsICJsYXRpdHVkZSIpXSwgcGNoID0gMTYsIGNleCA9IDAuNSwgY29sID0gInJlZCIpCmBgYAoK